Chirality and Equilibrium Biopolymer Bundles 
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We use continuum theory to show that chirahty is a key thermodynamic control parameter for 
the aggregation of biopolymers: chirahty produces a stable disperse phase of hexagonal bundles 
under moderately poor solvent conditions, as has been observed in in-vitro studies of F-actin [O. 
Pelletier et al, Phys. Rev. Lett. 91, 148102 (2003)]. The large characteristic radius of these chiral 
bundles is not determined by a mysterious long-range molecular interaction but by in-plane shear 
elastic stresses generated by the interplay between a chiral torque and an unusual, but universal, 
non-linear gauge term in the strain tensor of ordered chains that is imposed by rotational invariance. 

PACS numbers: 87.16.Ka, 64.75.-|-g, 61.30.Dk 



The structural integrity of cells depends on long, semi- 
flexible biopolymers-such as filamentous actin (F-actin), 
microtubules, and intermediate filaments-that make up 
the cytoskeleton. The controlled assembly and disas- 
sembly of these polymers is central to the division, de- 
velopment, and locomotion of cells [1]. In-vitro studies 
of their thermodynamic, structural, and elastic proper- 
ties have contributed importantly to our understanding 
of the physical mechanisms that underlie the function- 
ing of the cytoskeleton. For example, F-actin filaments 
adopt characteristic structural motifs inside cells: open, 
branched networks and filament bundles Q. Growing 
networks play an important role during cell spreading 
while dense bundles are a necessary component of cel- 
lular force generation (stress fibers) and cell protrusions 
(e.g. microvilli and stereocilia) . Networks and bundles 
appear as separate regions in the phase diagram of solu- 
tions of actin proteins and linkers under in-vitro condi- 
tions, [3, 0, 0, @| and much is now understood about the 
function of specialized linker proteins in the structural 
organization of actin filaments. Interestingly, even non- 
structural biopolymers like DNA 0] can exhibit these 
morphologies in solutions with generic linker units, such 
as polyvalent ions. 

The observation of stable biopolymer bundle aggre- 
gates poses an interesting puzzle. Equilibrium solutions 
of (charged) rods of fixed length are pi'edicted to form 
networks at low linker concentrations [8| . Beyond a crit- 
ical linker concentration, cylindrical aggregates should 
form without a limit on their diameter: bundles can 
grow in size until the supply of free rods is depleted. 
Explanations for the observed finite size of the bundles 
have focused on the possibility of long-range interac- 
tions [9], on kinetic limitations to the bundle size [lo| . 
and on assembly defects of toroidal-shaped biopolymer 
aggregates [lH. A recent study [l3| of fiuorescently la- 
beled F-actin mixtures has cast doubt on these expla- 
nations: non-toroidal bundles were found to maintain a 
well-defined, fixed diameter~in the range of 40-100 nm- 
that was greatly in excess of the electrostatic screening 
length even when filament exchange between and along 



different bundles was rapid. 

It is the aim of this Letter to argue that chirality is 
a key thermodynamic control parameter for biopolymer 
bundle aggregation. Chirality already is well known to 
play a determining role for the phase diagram of solu- 
ble chiral biopolymers. It is, for example, responsible for 
both the cholesteric liqrud crystal and the "blue" phases 
of concentrated DNA solutions [l3|- Here, we show that 
chirality also plays a key role for biopolymers in poor 
solvent conditions: in the presence of chirality, a stable 
disperse phase of twisted, finite-sized bundles appears, 
separating the bulk aggregation phase from the single 
filament (or network) phase. The large characteristic di- 
mensions of these chiral bundles is not determined by a 
mysterious long-range molecular interaction but, instead, 
by in-plane shear elastic stresses generated by the com- 
petition between a chiral torque and a non-linear gauge 
term appearing in the strain tensor for ordered chains. 

To demonstrate these claims, we will not construct an 




FIG. 1: Top view (left) and side view (right) of twisted bundle 
of hexagonally-ordered filaments. For a torsional strain, fl , 
the outer most filaments make an angle 9 = QR with the 
central axis. Hence, elastic strainchain bending and lattice 
shearbuild up as twisted bundle grows. 
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explicit model of F-actin based on specific assumptions 
about molecular structure, instead, we will rely on sym- 
metry arguments to construct the generic free energy of 
a hexagonal bundle of chiral, semi-flexible polymers of 
fixed length, L, where chains can slide freely along the 
local hexagonal axis. Deformations of the array are de- 
scribed by the in-plane chain displacement field u_l (r^ , z) 
of chains located at r_L (in the plain) as a function of 
height, z. Local chain-tangent unit vectors are given by 
t ~ 1 -f dzU± and local curvature of the chains is given 
by K = |9^u_l|. The rotationally-invarient in-plane strain 
tensor of the array is determined by uj^ through 

= ^ {diUj^j + dju^i - d,uj_ ■ (9jU^ - tit^) , (1) 

where indices refer only to in-plane directions. The first 
three terms are familiar from the classical theory of the 
elasticity of two-dimensional solids. The last term is a 
gauge term imposed by the chain connectivity: a uni- 
form rotation of the chains away from the bundle axis 
reduces their mutual perpendicular separation, generat- 
ing strain ^14] . The elastic strains captured by the gauge 
contribution to Uij play a special role in the following. 

The deformation energy is the sum of an in-plane en- 
ergy Ti.± that depends only on this strain tensor. For 
a chiral bundle with hexagonal symmetry, TC± can be 
expanded to second order in Uij as: 



(2) 



plus an out-of-plane term Tit that can be expanded to 
second order in t as 
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+ K2{Vi_ X t - z/Pf + K:i{d,tf^ . (3) 

The in-plane and out-of-planc contributions to the de- 
formation energy are linked, both by the relation t ~ 
z + dzVLi_ between the tangent vector and the in-plane 
displacement, and by the non-linear gauge term that ap- 
pears in the strain tensor eq. ([1]). In the conclusion we 
will discuss estimates for the Lame constants and /i^ 
and the Frank constants ii'i.2,3- In what follows we as- 
sume, ^ Ki^2 ■ Below we focus on bundles of variable 
radius R and fixed length L, such that L ^ R. 

The term K2P~^{'S/± x t)^ has the form of an exter- 
nal torque that tends to twist a bundle. We will show 
elsewhere that the resulting in-plane displacement field 
can be represented, to a very good approximation, by 
a pure torsional deformation where each cross-section of 
the bundle is rotated by rigid-body rotation over a cer- 
tain angle with respect to the sections above and below 
it (see Fig. [1]). Near the cross-section at z = , such a 
torsional displacement field can be written as: 



r(zfl)^ 

Ux(»', z) ~ r(zil)0 — f 



to second order in the torsional strain, fl . The defor- 
mation energy per unit volume E^iast (R, ^) for given 
and bundle radius R can be obtained by inserting eq. (j4]) 
into Ti± and Tit : 



E,iast{R,n) = fi^inRY 



. (5) 



Here, flj_ = 25Aij_/216, K2 = 2K2, K3 and 7 = 2K2/P 
(we have assumed the array to be incompressible). We 
assume the magnitude 9 = fii? of the strain to be small 
compared to one, as we are neglecting higher-order terms 
in Ti,±. The last three terms of eq. ([5]) result from a 
straightforward evaluation of ([3]) while the first term, 
which is non-linear and dominates for large R, is the 
shear stress produced by the gauge term in eq. ([2]) . Fig- 
ure [1] shows the physical origin of this term. Chains far 
from the central axis are increasingly tilted away from the 
bundle axis by the torsional strain field, leading to com- 
pression of the local hexagonal packing of chains along 
direction, f x t. The term, — 7^^, in eq. ([5]) acts as a 
chiral contribution to the aggregation energy density of 
the rod, with the first three terms the unavoidable energy 
penalty associated with chiral aggregation. There is, of 
course, a conventional achiral contribution as well to the 
aggregation energy density, which we denote as — e. 

For K2 <C K3 , the elastic energy is minimized by a tor- 
sional strain that depends on the bundle size as ^o{R) = 
ij/Afl^y/^R-y^{R^ + XlJ-'/^ where Agx = v^^^^ 
is an important characteristic length. The total free en- 
ergy of the bundle is the sum of the elastic energy eq. 
([5|)-with n = ^lo{R) -the achiral aggregation energy and 
a surface energy associated with the reduced binding en- 
ergy of filaments located on the surface of the bundle 
exposed to solution: 



AF{R) 



= -fo 



i?2)l/3 



nR^e + 27ri?S . (6) 



Here, fo = 37r(7^/16/2j_)^/'^ and E is the surface energy 
per unit area. 

Under conditions of thermodynamic equilibrium, the 
fraction P{R) of chains that will aggregate into bundles 
of radius _R-the cluster size distribution-in chemical equi- 
librium with the surrounding solution follows the Boltz- 
mann distribution llTII: 



P{R) oc exp 



AF{R) - hpqttR^ 



(7) 



(4) 



Here, fj, the chemical potential of filaments and po the 
area density of chains. In the L — > cx) limit, this cluster 
size distribution is dominated by the minima of AF{R). 
Near the point of bulk condensation-when p ^ pc = 
—eL/pQ-a,nd for large surface energies, the only minimum 
is an i? = boundary minimum for which AF — 0, 
corresponding to individual chains in solution or small 
clusters with an exponential size distribution. As the 
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FIG. 2: The dependence of the free-energy of a chiral bundle 
on bundle radius, R, for the columnar phase of aggregates at 
the bulk condensation point {fi — fic)- When surface tension 
is greater than Ec the free-energy minimum occurs for _R = 
(single filaments). For E < Ec the minimum occurs at a finite 
radius, _Ro > ^3± ■ The inset shows the phase diagram for the 
columnar phase of aggregates as function of filament chemical 
potential and surface tension. 



surface energy is reduced, a second minimum develops 
at i? = i?o(S) with a width SR cx L^^/^. For large L, 
an appreciable fraction of the chains can only condense 
into the clusters associated with this minimum provided 
that AF{Ro) — fiponRQ < 0, otherwise the boundary 
minimum at i? = dominates the distribution. Hence, 
a thcrmodynamically stable disperse phase occurs when 
the surface energy drops below a critical value Sc. When 
/i = /Zc, Ec is given by: 



E, = 



/o 



(8) 



For S < Ec , AF(i?o) is negative and the second mini- 
mum entirely dominates the cluster distribution function 
in the large L limit, and this corresponds to the stable 
dispersed phase of bundles (see Fig. [J). The peak po- 
sition shifts to values larger than Asj^ as the surface en- 
ergy is decreased, and in the limit of low surface energies 
it grows as Rq ~ ^'^/jl±Y,^. Thus, Asj, corresponds to 
a minimum bundle size at the bulk condensation point. 
Note that this minimum bundle size is determined by 
the non-linear contribution to the elastic energy, and it 
diverges in the fl± — > limit. Away from the bulk con- 
densation point, for /i < /ic, the minimmn in AF(R) — 
fipQirR'^ is reduced by the bulk free-energy penalty, lead- 
ing to reduction of the critical surface tension. Suffi- 
ciently far from the bulk condensation point we find a 

critical surface tension, Ec ~ /q^^3~^/^A;^jJ;(^c — 

The phase behavior of columnar aggregates is shown in 

Figured (inset). 

If the aggregate adopts full three-dimensional posi- 
tional order, then chains cannot slide past each other 



and a new term determined by the out-of-plane shear 
modulus must be added deformation energy density: 

EeiastiR, r!) = fi\\{nRf+fi^{nR)''+K3n^R^-jn , (9) 

where = Ml|/4 a-nd we have here dropped the term 
proportional to K2. This new term is actually familiar 
from the classical theory of elasticity. It is the strain en- 
ergy density of an elastic rod under uniform torsion [3| ■ 
In the physical range, where the magnitude of the elastic 
strain 9 — flR must be less than one, the new harmonic 
term fj,\\0'^ dominates-for /2|| « /2^-over the anharmonic 
term p,±9^ that previously played such a central role. The 
form of P{R) for a bundle in the solid phase is somewhat 
more complex, but the key features can be illustrated by 
analyzing the limit of large R. 

If we minimize only the harmonic (first) and chiral 
(last) terms of eq. ([9]) with respect to ft, then the opti- 
mal strain and the minimum energy density depend on 
R as do{R) ~ 7/2/2||i? and Eeiast — -"f'^ / 2fi\\R^ , respec- 
tively, where we must demand R S> 7/m|| in order for 
6* <C 1. The two anharmonic terms, Kj,B^ jR? and fl±9'^, 
are then both small compared to — 7^/2/i||i?^ for large R. 
However, in this regime the lowest order term contributes 
only the negative constant, — 7r7^/2/i||, to AF{R). While 
this constant is important in terms of competition of chi- 
ral aggregates with individual filaments, it alone leads to 
no preferred bundle size. Thus, to harmonic order, elas- 
ticity theory of three-dimensional chiral polymer solids 
predicts no maximum in P{R) . However, if we include 
anharmonic corrections, for sufficiently low surface ener- 
gies and near to bulk condensation (fj. > /Uc), then there 
is a negative minimum AF(R), and thus a dominant 
maximum in P{R) near Rq ^ {K3/T,y/^{-f / p,^^)^/^ when 
Ro < X3± and Rq - (m±/E)1/3(7//I||)4/3 for i?o > Ag^ . 
Importantly, in both hexagonal-columnar and hexagonal- 
solid phases, preferred bundle sizes are only possible if the 
effects of anharmonic elasticity are taken into account. 

To be clear, our approach is generic, relying on symme- 
try arguments to dictate the form of the elastic energy of 
ordered chains, rather than a detailed, microscopic de- 
scription of intermolecular forces between biopolymers. 
Nevertheless, in order to apply the continuum theory to 
F-actin bundle formation, we need estimates for the elas- 
tic moduli and the Frank constants for realistic biopoly- 
mer bundles in solution. For a bundle held together 
by linker units-proteins, cations, ete.-the Lame con- 
stants A_L and p,±_ can be crudely estimated to be of 
order AE/lcP , with A£: - 1 - lOfcsT the binding 
energy per linker unit, £ the spacing between linkers 
along the polymer backbone, and d ^ 1 — 10 nm a 
typical molecular length (e.g., the spacing between the 
chains). The respective Frank constants Ki and K2 for 
splay and twist will be estimated as AE/£ . The Frank 
constant for bend is determined not primarily by link- 
ers, but by the bending stiffness of the individual chains 
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through through Kj, ~ ksTEp/d? with £p the persistence 
length. For F-actin tp is of the order of 10 /xm so in- 
deed i^s ^ Ki 2 for F-actin filaments as claimed earlier. 
Given these estimates, the minimum size of bundles in 
the hexagonal columnar phase at the onset of aggrega- 
tion, As^ ^ J dl 



/dip ~ 100 nm, indeed is of the order of 
the observed bundle size [13] ■ The surface energy will 
be estimated by E ~ Ai?*/^d, where, in general, Ai?* 
can be less than AS, for example, when surfactants are 
absorbed on the bundle surface. The ratio Ec(//c)/E is 
of order {d/ P)'^/^{d/ X3±y/^ AE/ AE* . When the dis- 
perse columnar phase is stable, this parameter must be be 
larger than one. Disperse bundles are thus stabilized by 
reduced persistence lengths, by the action of surfactants, 
and by a short preferred chiral pitch P, preferably on the 
order of molecular scales. For large persistence lengths, 
the radius of disperse bundles with solid order scales as 
Rq - (A'g/E) 1/5(^/^11)4/5 _ ^1/5^ ^j^j^j^ considerably 
reduced in comparison to the scaling of minimum size for 

1/2 

columnar bundles, Asj, ~ £p . This indicates the radii 
of solid bundles is much less than that of the columnar 
phase, owing to the increased resistance to torsion [l9| . 

Confirmation of the description presented in this paper 
demands the observation of torsion in bundles of chiral 
biopolymcrs by microscopy. Cryo- AFM images of F-actin 
filaments by Shao et al. provided the first clear evidence 
of twisting More recently, Ikawa et al. carried out 
a careful AFM study of F-actin filaments for different 
Mg^+ ion concentrations [2l[. For low Mg^+ ion con- 
centration, a network is observed which breaks up into 
braided bundles of a small number of filaments at in- 
creased linker concentrations. Finally, synchrotron X-ray 
studies of F-actin bundles show small but clear devia- 
tions from perfect hexagonal packing 0| , which would 
be consistent with the anharmonic chiral strain discussed 
in this paper. The combination of these observations and 
the theoretical analysis of the present paper we believe 
provides mounting evidence for chirality as a determin- 
ing factor for the formation of thermodynamically stable 
F-actin bundles of finite size. 

We conclude by noting that the present analysis did 
not allow for the possibility of topological defects. Kamien 
and Nelson 2^ (KN) showed that an infinite hexago- 
nal columnar array of strongly chiral polymers is unsta- 
ble against break-up into finite-sized regions of localized 
double- twist. In this so-called "moire phase", double- 
twisting bundles are arranged in a triangular lattice, with 
neighboring bundles meeting along a tilt-grain boundary, 
necessary to accommodate the mismatch in chain orien- 
tation between adjacent bundles. The size scale £,m of 
the bundle is determined by the repulsive energy between 
screw dislocations along the grain boundary and the in- 
trinsic torsional stress as discussed in the present paper. 
We speculate that the KN moire phase could be viewed 
as the result of bundle-bundle aggregation produced at 



low surface energy. As E is decreased, and the equilib- 
rium bundle size Rq approaches from below, individ- 
ual merge, hiding surface exposed to unfavorable solvent 
contact by introducing localized grain boundaries where 
neighboring bundles meet. Alternatively, we may view 
the moire phase as a precursor to the chirality-induced 
break up of the hexagonal columnar phase. In a future 
study, we aim to extend our analysis to the more general 
case including screw dislocations as well as to the case of 
the toroidal geometry observed for bundles of long DNA. 

It is a pleasure to acknowledge many useful discussions 
with R. Kamien, S. Safinya and M. Henle. This work was 
supported by the NSF under DMR Grant 04-04507. 
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